Coherent electron tunnelling dynamics in polyacetylene 

D. Psiachos 

Abstract 

We examine, theoretically, using a mixed classical-quantum electron-lattice dynamics model, 
electron transfer in a molecular junction system with a polyacetylene wire at its centre. We identify 
two qualitatively-different transfer regimes: hopping and tunnelling. We discuss the criteria for 
achieving each one and for minimizing inelastic scattering and decoherence arising from the coupling 
to the lattice, and we connect our main results to quantities derived from electron dynamics in 
simpler, three-state model systems. We identify the requirements to have near-ballistic transfer. 

PACS numbers: 73.22.-f,73.63.NM,63.20.kd 
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I. INTRODUCTION 



Electron transport through organic molecules or chains of molecules has been an active 
area of fundamental research because of its potential applications to molecular electronics. 
Recently, measurements of the conductance of single molecules or atomic wires attached to 
leads have become very refined^. A simple expression for the electronic conductance through 
a molecule is provided by the Landauer formula, which states that the conductance between 
two leads is proportional to the transmission, which is assumed to be constant for the narrow 
energy range covered by the leads. The actual behaviour of the transmission depends on the 
scattering processes involved. Some of the methods used to describe the scattering include 
Green's function methods, transfer matrix approaches such as the electron scattering quan- 
tum chemistry (ESQC) method^, or the use of Lippmann-Schwinger scattering operators 
for mapping the non-equilibrium system onto an equilibrium one^. All of these methods are 
time- independent, thus enabling the use of the Landauer or, in the case of several leads, the 
Landauer-Biittiker formula, to evaluate the current. Time-dependent methods include non- 
equilibrium Green's function-based approaches or approaches based on the time- dependent 
Schrodinger equation applied to the key parts of the system, i.e. the tunnelling electron, 
and enable a direct evaluation of the transmission and the current. 

The treatment of the electron-lattice interaction within semi-empirical formalism assumes 
the validity of one of two regimes: strongly or weakly-localised electron-lattice interaction. 
The Marcus theory^ treats electron transfer non-adiabatically with a rate expression in- 
volving the initial, donor, and final, acceptor, states and geometries expressed in terms of 
two key quantities, presumed to be separable owing to the slow intermolecular vibrations 
compared to the fast electron transfer rate: the difference in energy of the ground-state con- 
figurations before and after electron transfer, and the non-adiabatic contribution, called the 
reorganization energy, which is defined as the electronic excitation energy of the initial state 
in its ground-state configuration. As in the Marcus assumption of distinct donors/acceptors, 
the Holstein model^ also treats the electron as being highly localised in space. These tech- 
niques are most appropriate for non-homogeneous systems or systems containing defects or 
other strong scattering centres where there is a clear distinction between donor and acceptor 
and the transfer is accomplished by hopping between these two species. In more homoge- 
neous systems, there tends to be a more delocalised electron-lattice interaction, leading to 
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relatively adiabatic electron transfer. 

The different models for electron transfer have been developed with the goal of describing 
polaron motion in oligomers. The Holstein model describes "small-polaron" formation, 
or localised interaction with a particular molecule (site) and the electron conduction is 
accomplished by distinct hopping, a non-adiabatic process. The structural perturbation 
is similarly highly localised. Another regime, the "large-polaron" model, which leads to a 
longer-ranged structural distortion and at the same time, a long-ranged electronic effect, 
is suitable for periodic systems. The best-known description for large-polaron formation is 
the SSH^~— model. Polaron transport in oligomer chains using the SSH model driven by an 
electric field has been studied extensively^— . The main result of these studies is that at 
high field strengths polaron dissociation occurs because of a complete decoupling of electron 
and lattice motion time scales for highly-energetic electrons, which leads to increased charge 
mobility as the carrier becomes a free electron. 

Recently, the Holstein and SSH models have been combined, in order to describe both 
the inter- and intra-molecular interactions respectively of a crystal formed of pentacene 
chains^. While the inter-molecular interactions, originating from Van der Waals forces, are 
weak, they are not neglected because in realistic molecular crystals there is a high degree 
of disorder in the molecules' arrangements, and this was found to play a significant role in 
determining the type of charge transport at nonzero temperatures - going from bandlike to 
diffusive hopping transport as the disorder increases. 

In our work, we work with a single chain of atoms only, and as such within the "large- 
polaron" regime. Specifically, we study the influence of electron-lattice coupling on electron 
transfer in a polyacetylene (PA) chain. Rather than using an electric field, we consider 
electrons coming from a reservoir of varying chemical potential, and study the influence of 
this as well as the coupling to the PA chain on the electron-lattice dynamics timescales and 
the charge propagation through the chain. 
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II. METHODOLOGY 



A. The SSH Hamiltonian 

We use the SSH Hamiltonian^ to describe the electronic and structural properties of the 
PA chain. This model uses a basis of site-centred ir orbitals in the tight-binding approxima- 
tion to describe the electronic properties, while the lattice is treated classically. Regardless 
of the number of electrons, the basis size is fixed to the number of atomic sites, as per the 
LCAO approximation. Explicitly, the Hamiltonian is split into electronic and ionic parts: 

H = H e i + H ion (1) 

where 

H e l = — ^ ] tn+l,n ^ C n+l,o- c n,<r + C n,a C n+^-,'^) + ^ ] ^n^ n .a C n,rj (2) 
n,cr n,cr 

is the electronic part where, the operator c n a adds an electron of spin o to site n and where 
in the SSH model, the site energies e n are set to zero and 

t n +i,n = t - a (u n+1 - u n ) (3) 

describes the hopping of the electron. The hopping is modified by the electron-phonon 
coupling parameter a according to the displacements u n of the ions, from their positions in 
a perfect, undimerised chain and its range is confined to nearest neighbours. 
The lattice part, 
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H ion = V ion + f ion = - ^ K (Un+l ~ U n f + ^ Tfh ( 4 ) 



consists of a harmonic potential energy operator with spring constant K and a kinetic energy 
operator. The set of parameters used is standard in the SSH model literature as applied 
to PA chains^: t =2.5 eV, a=4.1 eV/A, K=21 eV/A 2 , and M=1349.14 eV fs 2 /A 2 . The 
electronic and ionic Hamiltonians are thus linked by the ionic displacements u n . 

We include the site-energy term in Eq. |2] because we will consider a chain where each end 
is coupled to a molecular level with energy e en d- All other e n are set to zero. In addition, 
we will modify the coupling between the molecular level and the chain to a value t en( i. We 
will take the molecular levels and couplings at both ends to be the same. A schematic of a 
chain with its end sites coupled to molecular levels is shown in Fig. [U The coupling to the 
electrode t en d is fixed and does not depend on the distortion of the chain's end sites. 




FIG. 1: Chain with nearest-neighbour hoppings based on the site distortions u (see Eq. [3]), coupled 
to molecular levels. The total number of sites including the molecules at the ends is N. 

B. Solving the coupled ion-electron system 

We use the Ehrenfest approximation^^—, a non-adiabatic semi-classical scheme 
in which the ions undergo classical motion while the electrons are treated quantum- 
mechanically, evolving in time. The adiabatic approximation, also known as Born- 
Oppenheimer (BO) dynamics, a simpler scheme in which the electron orbitals are treated 
time-independently, affected only by changes in the ionic potential, and thus not deviat- 
ing from their initial superpositions of Hamiltonian eigenstates, would not lead to electron 
transfer which is what we want to study here. To that end, the equation of motion of the 
electrons \i/)k(t)) of state k is given by 

ih djm_ = Aa m)) (5) 

while that for the nth ion is 

k,occ 

which we integrate using the velocity- Verlet algorithm 

r 2 

u n {t + t) = —F n + 2u n (t) - u n {t - t) (7) 

for a sufficiently-small timestep r with which to maintain energy conservation - e.g. a r 
of 0.02 fs conserves the total energy to within 0.13 meV/site over 1 ps. In Eqn. [6], we use 
the definition of the Hellman-Feynman force because it is this and not the derivative of the 
negative of the electronic energy which yields the expectation values of the ionic trajectories. 
In practice however, we found virtually no difference between the results obtained by the 
two force expressions for systems as large as the ones we studied. 

Letting the electronic eigenstates arising from a solution to the time-independent 
Schrodinger equation 

H e l \Xm) = e m \Xm) , (8) 
5 
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where H e \ implicitly depends on the ionic configuration {u n }, be represented by the set {Xm}, 
the wavefunctions of all the electrons k can be written as superpositions of the eigenstates 
as 

^kit)) = S ^g m A t )\Xm) (9) 

m 

The initial conditions of the problem are divided conceptually into two types: the first 
type of initial condition is g m ,k{Q) = $mk and is imposed on the N wavefunctions of the N-site 
system, while the extra electron is set up in a superposition state in general. Thereafter, 
all of the ipk{t) evolve according to Eq. along with the boundary condition that they are 
continuous across time-steps, and become non-stationary states, written in terms of the 
eigenstates corresponding to the current {u n }. 

In all of the results presented in this paper, we treat all of the electronic wavefunctions 
with the Ehrenfest method. Before an extra electron is added to the system, the ground 
state of the neutral chain, whether coupled to a junction or not, is determined. This is done 
via the procedure described above with the ionic velocities held to zero. Here, the ends, 
which may correspond to coupled molecules, are free to move while a total length constraint 
is imposed on the system. The purpose of this is to obtain displacement values for the ends 
which are consistent with the parameters of the Hamiltonian, thereby minimising the finite- 
length effects. Afterwards, for the dynamics, the length constraint is removed and the end 
states are held fixed to their relaxed values because keeping the length constrained instead 
of fixing the ends led to the plots of displacement u versus site number being tilted, as the 
chain would translate, and this made the results difficult to compare and analyze. Imposing 
no boundary conditions whatsoever results in highly-distorted ground-state structures due 
to the finite length of the chains, obscuring the phenomena which we would like to study. 

III. SYSTEMS STUDIED 
A. Isolated PA chain 

Two types of open neutral chains were studied as starting points: with an even or an 
odd number of sites. The ground state of the even chain leads to the well-known dimeri- 
sation pattern^, shown in Fig. Hk, while that of the odd chain contains a soliton, again as 
expected^. 
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After obtaining these ground states, an electron was prepared in either the LUMO, or 
in a non-stationary state at one end, and the dynamics were performed. In the even-site 
chain, the case of the LUMO preparation leads to a polaron^^^, a structural deformation 
consisting of a soliton/anti-soliton pair moving in tandem with the electron localisation. We 
focused on even-site chains in order to avoid considering soliton dynamics along with the 
transfer of the extra electron. The effect of the Ehrenfest approximation as opposed to the 
BO one is clear from Fig[2j The energy of the extra wavepacket, initially prepared in the 
LUMO, closely follows the LUMO but nonetheless deviates from it at times. 
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FIG. 2: The energy of the added electron deviates from the LUMO, in Ehrenfest dynamics. 

For the electron prepared at one end, the charge density oscillated back and forth (see 
Fig. |3h). This oscillation was virtually all due to the extra electron - the other electrons in 
the chain had practically no contribution to it. However, the chain distortions were very 
large, several times greater than those corresponding to the dimerisation (~ ±0.04A, see 
Fig. UK), and in practice, would require anharmonic potentials to treat, and may lead to 
chain breakup. It is important to note that as the electron wavefunction oscillates back 
and forth, it distinctly passes through all sites between the two ends. Therefore, this is a 
hopping transfer regime as opposed to a tunnelling one. 

In order to determine the effect of the lattice on the electron propagation, we compare 
the propagation with that for the case of a frozen lattice (Fig. [3b). In a frozen lattice, 
the propagation is quasi-periodic, owing to the non-commensur ability of the eigenvalues in 
the eigenstate superposition comprising the wavefunction of the propagating electron. For 
increasingly longer times, according to the Poincare theorem applied to quantum superpo- 
sitions of eigenstates with non-commensurate eigenvalues, the oscillations will approach the 
original value increasingly faithfully and eventually lead to a value near 1. This is the reason 
why the value of 1 is not reached during the few ps of the simulations for the frozen ion case. 
Beyond that, the electron-lattice interaction destroys the coherences (off-diagonal entries of 
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the density matrix), which manifests itself in a decrease in site populations in Fig. |3]3 for 
the simulation with full dynamics, as compared with the frozen-ion case. The off-diagonal 
entries of the density matrix appear in the electron density due to the extra electron as a 
result of intereferences from the eigenstate superposition. 
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FIG. 3: (a) The populations of the two end sites 1, 60 and the middle site, 30, for electron-ion 
dynamics for a chain not coupled to molecular junctions where an extra electron was added at site 
1 at t=0. (b) Comparison with the frozen ion case. The population shown is only due to the extra 
electron (see text for an explanation). 



B. Chain within molecular junctions 

The molecular junction coupling parameters t en d and e en d were varied across a wide range 
of values in order to search for an optimal set which would lead to low electron decoherence, 
meaning a small influence of the lattice on the electron transfer. This required first obtaining 
a relaxed ground-state structure for the lattice, and then performing dynamics using the 
added wavepacket. We only performed dynamics for cases which looked promising based 
on the frozen-lattice transfer calculations. The ground-state structure for a chain coupled 
to a molecule at each end deviates significantly from that of a dimerised chain (Fig. Hh)- 
Most notable, apart from the end effects, are the kinks introduced at large e en d and small 
tend, (Fig. Ht) which migrate towards the ends for smaller e en d and larger t end (Fig. HJd). 
Usually, the site population over time can be approximately represented by a single Rabi 
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FIG. 4: (a) Ionic distortions u(n) in Afor a chain with no coupling to a molecular junction at its 
ends, (b) e en d=0.3 eV and t en d=l-0 eV, and (c) e en d=0-8 eV and i en d=0.4 eV. 

frequency, the frequency of oscillation for an electron in a three-state quantum system, by 
a projection technique^. Characteristics such as short Rabi oscillation timescales with a 
single dominant Rabi frequency, usually led to the most unaffected by the lattice vibrations 
transfer (see later in Sec. II II CI) . The reverse is found to be the well: the electron 

oscillation did not affect the structure a lot, which is a different result than the one found 
by Ness and Fisher— who used a model with quantized phonon modes. What completely 
destroyed the electron oscillations was the presence of noise which competed with the main 
Rabi frequency and led to destructive interference in the oscillations. 

Regarding the charge density oscillations, as with the case of the chain not coupled to any 
junction, only the extra electron reached the other side. The other electrons nevertheless 
caused a significant amount of charge- density oscillation, which was correlated with the 
evolution of the lattice deformation. In Fig. [5] we show a case also in the "hopping" regime, 
albeit not as extreme as the case of Fig. [3] which had equal occupation of all sites. Still, there 
is a non- negligible occupation of all of the interior chain sites (only a central site is shown 
in the figures), and as discussed below, numerous Rabi frequencies, which do not show up 
at the small times in which coherent oscillation is maintained. 

It is important to note that distinct hopping across the chain by the extra electron was 
not always found to occur - we also identified another limit: that of tunnelling. In Fig. [6] we 
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FIG. 5: (a) The populations of the two end sites 1, 60 and a central site, 30, for electron-ion 
dynamics for a chain coupled to end molecules with energy e en d=2.0 eV and coupling to chain 
t en d=l-l eV where an extra electron was added at site 1 at t=0. (b) Comparison with the frozen 
ion case. Only the contribution from the extra electron is shown (see text for an explanation). 

show a case in the "tunnelling" regime which had relatively low scattering of the electron by 
the lattice. It is called tunnelling because there is almost no occupation of the chain sites, 
with one dominant electron oscillation frequency connecting the defect levels arising from 
the end states (see Sec. IIII CI) . 

An interesting case, where the electron oscillation frequency changes dramatically with 
the dynamics, is shown in Fig. [7J Also, the oscillations are incredibly persistent over very 
long times. Generally, the electron oscillations are stable over long times throughout the 
tunnelling regime because the electron tunnels between the two fixed end molecules with 
minimal impact from any structural distortions which may be occurring in the chain. 

C. Rabi frequencies 

We systematically studied the oscillation of an electron in a frozen lattice by varying 
different fundamental quantities: the lattice vibration parameter K, the electron-phonon 
coupling parameter a, and the electronic parameters t en d and e en d of the end molecules. We 
determined the contribution of the pair of eigenstates comprising each Rabi frequency to 
the electron oscillation by projecting the eigenstates onto the end-molecule wavefunctions^. 
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FIG. 6: (a) The populations of the two end sites 1, 60 and a central site, 30, for electron-ion 
dynamics for a chain coupled to end molecules with energy e en d=0-6 eV and coupling t en ^=0.4 eV 
where an extra electron was added at site 1 at t=0. (b) Comparison with the frozen ion case. Only 
the contribution from the extra electron is shown (see text for an explanation). 

Explicitly, this weight for the (n,m) pair of eigenstates is equal to 



D 



nm — Wl 



P, 



nd) ((fiend An 4>1 ) (10) 



where the operator P, = \Xi)(Xi\ is the projector of eigenstate i and the two electrode sites 
are labelled (fii and <fi en d- The quantity in Eq. [TO] is actually the maximum value, as the 
weight oscillates in time according to the Rabi frequency. In three-state quantum systems 
there is a relationship of inverse proportionality between the Rabi frequency and the weight 
(see Eq. 3 in Ref. |22| ) which we exploit in our analysis of Sec. IIVBI We focused on varying 
the end-molecule parameters t en d and e en d in our analysis of the Rabi frequencies. 

The hopping case (Fig. [3] and [5]) is characterised by the presence of a large number of 
competing Rabi frequencies while in tunnelling (e.g. Fig. [6]), just one or two frequencies 
dominate the oscillations. However, even in hopping, at small times there seems to be 
just one main frequency present, corresponding to a fast time scale, on which there is 
superimposed noise from an even faster one. The other frequencies, with a longer period, 
eventually show up, and lead to an indistinct oscillation but by that time, the oscillation is 
already destroyed by the electron-lattice interaction. In Sec. IIV Al we discuss which choices 
of parameters lead to these different regimes and how the transmission occurs. 
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FIG. 7: (a) The populations of the two end sites 1, 60 and a central site, 30, for electron-ion 
dynamics for a chain coupled to end molecules with energy e en ^=0.7 eV and coupling t en ^=0.4 eV 
where an extra electron was added at site 1 at t=0. (b) Comparison with the frozen ion case. Only 
the contribution from the extra electron is shown (see text for an explanation). 

IV. ANALYSIS 
A. Defect states 

The perfect dimerized chain has a double bond at both ends and an energy gap of 1.39 
eV. When t en d is significantly different from to, defect states appear in the gap. If the chain 
is still well-dimerized, as for the case of Fig. Ht, there are only two defect states, one filled 
and one empty, degenerate and at nearly the energy e en d] otherwise there are two additional 
defect states, split off from the valence and conduction band respectively, consistent with 
a bi-polaronic distortion such as that shown in Fig. Hfc. The latter case is the most usual, 
as the case of only two defect states is confined to marginal cases (very small e en d\ Fig. M 
shows the one-electron energy level structure of the system attached to end molecules, for 
the most common case, that of a bi-polaronic distortion. For low t en d the empty-state defect 
levels are nearly degenerate and approximately equal to e en d- The bi-polaronic distortion is 
gradually lost for t en d > e en d and progresses to the ends for t en d S> e en d, and for sufficiently 
large t en d/e en d, there are no defect states at all in the gap. 

The closeness in energy of the empty states for low t en d means that the Rabi oscillation 
period is very long. This is also the pure "tunnelling" regime because there are virtually no 
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FIG. 8: Schematic of the energy levels for a system with two kinks in the distortion u(n), e.g. 
Fig. 2b. In a case such as Fig. 0b, there is no longer a dominant Rabi frequency and transmission 
takes place through many levels, including those in the conduction band (CB), filled defect levels, 
and valence band (VB). 

other contributions from other pairs of eigenstates to the electron oscillation between the 
ends. 

For t en d *C e em i, the defect levels are present, owing to the bi-polaronic distortion through 
the lattice, but only the filled ones are in the gap while the empty ones are located in the 
CB, assuming that e en d is large enough that it is in the CB. A very low t en( i combined with 
a large e en d leads to a very small Rabi frequency, or very slow electron oscillation, which 
may be undesirable. Nevertheless, this is the pure tunnelling regime, as these defect energy 
levels are highly localised at the end molecules. However, as t en d increases, more CB states 
around the energy e en d suddenly become involved in the transfer leading to a noisy spectrum. 
Thus, a large e en d which leads to defect states in the CB is not very practical for achieving 
stable coherent oscillations. On the other hand, for t en d 3> e en d, there are no defect levels 
as the chain is dimerised and the transfer occurs through a large number of levels, both 
filled and empty. The extreme case is that of Fig. [3j For both of these situations, the 
system is in the hopping regime as the tunnelling through numerous energy levels implies 
that the wavefunction of the extra electron has a non-negligble projection onto the chain 
sites. However, for intermediate values of t en d and e en d, we are able to obtain tunnelling and 
a fairly large, single dominant oscillation frequency, which is stable under the effect of lattice 
dynamics for at least a few cycles, over a wide range of molecular-junction parameters. 
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B. Effect of dynamics on Rabi frequencies 



Here, we describe the effect of ion dynamics on the Rabi frequencies, in order to explain 
the decoherence mechanism. A recent study^ examined the decoherence mechanism for 
electron wavepackets added to chains and determined that for long chains, decoherence is 
caused by population transfer into other states due to the dense energy level structure. 
Our results, which show a decay in population of the initial wavefunction superposition, 
especially in the case of the defect levels located in the conduction band, agree with this 
explanation. 
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FIG. 9: Variation of key quantities with time for the case t en d=0A eV and e en d=0.7 eV. The top 
panel shows the eigenvalues, L: LUMO, LI: LUMO+1, L2: LUMO+2. The middle panel shows 
the inverse of the Rabi frequency, T(t) = l/f2(i) for the two competing Rabi frequencies formed 
by the two pairs of energy levels shown. The bottom panel shows the variation with time of the 
weights of these Rabi frequencies. For clarity, the data shown is a running average over 100 fs. 

In addition, we find that when multiple Rabi frequencies compete, they interfere. As 
well, in the hopping regime, the electronic oscillation is substantially degraded by the ionic 
dynamics, as mentioned earlier at the end of Sec. IIV Al However, for several cases of two 
competing Rabi frequencies, we find that the electron dynamics exhibit exceptional stabil- 
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ity, particularly when the two Rabi frequencies contribute roughly equally to the oscillation. 
The reason is that their weights oscillate out of phase with each other. In Fig. [9] we show 
the time dependence of the energy levels involved in the oscillation, the inverse Rabi fre- 
quencies derived from these energy levels, and the maximum weights of the latter, for a case 
which had great stability in the electron oscillation. As mentioned earlier in Sec. II II C\ the 
quantities T(t) and D(t) are correlated. We create a simple model for combining the two 
Rabi oscillations, omitting the time dependence from T(t) as it is already found in D(t). 
If we denote the average value of the ith inverse Rabi frequency (connecting two particular 
energy levels) as Tj = 1/fij, and its maximum weight Di(t) (see Eq. [10]), oscillating with a 
frequency uj , we may write 

= % 

D x (t) = Ih + Bxcos^ot) 
T 2 (t) = T 2 

D 2 (t) = D~ 2 - B 2 cos(uj t) (11) 

and the combined oscillation period becomes: 

^_ T 1 D 1 (t)+T 2 D 2 (t) 
Tnet{t) ~ D 1 (t)+D 2 (t) ■ (12) 

We found that this expression represents the electron oscillation period very well. The 
main cause of any deviation from the prediction of Eq. [12] and the data is due to other 
secondary inverse Rabi frequencies, which are smaller and tend to bring the actual value 
down. The time- dependence in Eq. [12] is generally quite weak, owing to the opposite signs 
in the time-dependent quantities of Eq. [TT] and this leads to the exceptional stability seen 
in Fig. [7] We determined the average value T av of T net (t) and compared it with the main 
oscillation period for several cases, which we show in Table [H Regarding the oscillation 
frequencies u , they mainly corresponded to the bare (without the extra electron) lowest 
phonon modes of the respective system but not always, and importantly, no relation was 
seen between the matching or not of the ojq with a bare phonon frequency and the stability 
of electron oscillations in the respective system. 
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TABLE I: Parameters of Eq. [TT]for several cases of e en d (in eV) for fixed £ en d=0.4 eV, as deduced 
from the dynamics, for comparison with the average value T av = T ne t(t) and with the main 
electron oscillation period T ac t ua i. Note that T(0) is the frozen inverse frequency. If two inverse 
Rabi frequencies compete, the values corresponding to both are given, the same one always in 
parentheses. All times are in fs. 

V. CONCLUSIONS 

In summary, we investigated the inelastic electron transfer in a polyacetylene wire within 
molecular junctions and we found that under certain conditions, which we identify, there 
is near-ballistic transfer, in two distinct regimes: tunnelling and hopping. The tunnelling 
case resulted when the defect levels from the electodes were in the band gap of the chain 
whereas when they are in the conduction or valence band, tunnelling occurs through multiple 
energy levels leading to a non-negligible occupation of the chain sites, a regime which we 
term hopping. Ness and Fisher— have calculated tunnelling electron transport for a chain 
coupled to metallic leads, using a model incorporating a quantized description of the main 
phonon modes, and find a clear indication of polaron transport, or a distortion of the lattice 
which propagates along with the electron. On the other hand, another study-, based on the 
SSH model with a classical, single-oscillator description of lattice vibrations, does not find 



1(3 



any polaron transport. It does find however, that energy levels arising from a pre-existing 
polaronic distortion aid in the tunnelling process. Our work suggests that not polaronic 
per se, but defect levels in the energy gap, produce tunnelling transfer. We expect that 
there will be dramatic changes in the transport characteristics, particularly for the cases in 
which the defect levels are located in the conduction band, if we consider electron-electron 
interactions, and this extension is presently under investigation. We do not expect that 
a full treatment of transport from the coupling of the system to metallic leads will lead 
to qualitative changes. While we find that multiple competing Rabi frequencies lead to 
hopping and a destruction of the electron oscillations because of their different behaviour 
under the dynamics, we find a large stable regime where two competing Rabi frequencies 
interfere destructively in order to yield exceptionally-stable electron oscillations at a faster 
rate than if the ions were completely static. Generally in the tunnelling regime, we predict 
that when connected to a source of current, the conductivity will increased as a result of the 
ion-electron interactions. 
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